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SECTION I 
INTRODUCTIO N 

The selection of the nozzle design for a hypersonic scramjet should be 
done based on a compromise of external and internal flow requirements, re- 
lated to drag, lift, and pitching moments of the vehicle, and thrust of the 
engine. in addition, structural and weight considerations should be included 
therefore, a complex design study is required based on engineering judgement 
that cannot be obtained solely from complex numerical programs. In the pre- 
liminary design stage, a logical sequence appears to be first obtaining a 
satisfactory aerodynamic design criteria, and then refining this by intro- 
ducing weight and structural considerations. The present work is directed 
toward developing criteria for selection and methods of analysis directed 
toward definition of a satisfactory aerodynamic design. 

The methodology employs three steps differing in levels of sophistication. 

STEP 1 - In this initial step, the effects of many parameters 
involved in performance are evaluated on the basis of 
a simplified wave analysis in order to determine their 
relative importance and their appropriate values. 

This first phase is devoted toward selecting the location, size and 
expansion ratio of the engine streamtube and the ratio between the expan- 
sion in the transverse direction to the expansion in the direction normal 
to the eng i ne. 

Figure (l) is a schematic of the important vehicle-engine parameters; where 
the lift is given by the suction force on the upper surface L^, by the spil- 
lage pressure , by the engine cowling pressure L^, (which is a function 
of Oy) and by the expans ion of the engine L^. The engine divergency in plan 
view o'^ induces a lift Lj-, and the divergency of the jet a ^ contributes to 
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FIGURE 1 


This step involves the use of a "source flow" charac- 
teristic analysis; approximate values of lift, thrust 
drag and pitching moments are obtained using the range 
of parametric values determined in Step 1. Additionally 
the influence of the nozzle on the external flow will be 
evaluated employing two dimensional and conical flow ap- 
proximat ions . 
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STEP 3 - As a result of Steps (l) and (2) the values of the main 
parameters are approximately known and preliminary per- 
formance levels have been obtained. Then, the most 
promising configurations can be analyzed in detail. 

A brief description of the simplified wave analysis of Step 1 is presented 
in Section II of this report. A numerical procedure incorporating the ideas 
outlined in this section has been developed and is described in detail in 
Reference (1). 

The "source flow" approach is briefly discussed in Section III. A complete 
description of this technique and its application to the analysis of scram- 
jet exhaust flow fields is contained in Reference (2), while a description 
of the numerical program and user manual is given in Reference (3). 

A description of the vehicle-engine interference is contained in Section IV. 
The analysis provides for determination of a preliminary exhaust plume shape 
based on quasi-two dimensional considerations as well as an assessment of 
the cowling lift and drag based on quas i -con ica 1 flow solutions. 

Steps (1) and (2) yield a preliminary estimate of the range of nozzle para- 
meters to be considered. After completion of this preliminary design work 
a full three dimensional calculation can be made using analysis and the 
program described in References (A), (5) and (6). This will permit further 
design refinements and the exact computation of thrust, drift and pitching 
moment for the vehicle. Thus, lengthly three dimensional calculations are 
avoided until simpler parametric or quasi two dimensional calculations have 
defined the basic nozzle configuration. The present approach allows the 
nozzle designer to simply investigate and assess the effects of various para- 
meters prior to employing three dimensional calculations. 


- 3 - 



TR 195 


SECTION I I 

PARAMETRIC APPROACH EMPLOYING SIMPLE WAVE TECHNIQUES 

The type of nozzle considered is depicted in Figure (2) and is seen to be 
a combination of several nozzles merging into a single nozzle. The in- 
dividual nozzles, before merging, expand the flow in two perpendicular 
planes as depicted in Figure (3). In this figure, the lateral expansion 
is depicted as occurring within two parallel walls. However, more general 
configurations can be considered which require a greater amount of lateral 
expansion and, hence, require diverging walls in the region of the indi- 
vidual nozzles as depicted in Figure (k) . 

In the preliminary analysis of Step 1, it is assumed that the jets, after 
merging, are bounded by sidewalls which extend downstream of the merged 
station. Then, each nozzle can be analyzed by considering two regions. 

The first region extends from the burner exit to the station at which the 
individual jets merge (Station A-A of Figure 3), while the second region 
is that downstream of the merged station. Note that the side expansion of 
the nozzles produces external compression. This lateral deviation is bene- 
ficial since for a given value of total expansion it decreases the deflec- 
tion angle of the exiting jet flow in the vicinity of the cowl. At the 
same time, the corresponding external compression does not yield too much 
of a drag penalty since the compression waves generated are reflected on 
the rear part of the wing resulting in lift and a decrease in drag due to 
the wi ng. 

A main consideration of this analysis will be to select the range of im- 
portant parameters to be considered. The initial flow (at the burner exit) 
is assumed known and is represented as an average uniform flow, obtained by 
properly averaging the actual flow at the burner exit station. The assess- 
ment'of changes introduced by flow nonuniformities at the entrance station 
is performed in Step 2 of the analysis. The parametric analysis is performed 
in conjunction with isentropic expansions from uniform initial conditions. 
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Initially assume that the sidewalls are straight and extend to the end of 
the nozzle, such that the flow is entirely contained within the sidewall. 
Additionally, assume that the internal vertical walls terminate along a 
line LN, parallel to the initial line AB which represents the burner exit 
station as depicted in Figure (5). 

It should be noted that the flow properties can be determined in separate 
steps due to the segmenting afforded by the waves of the flow field. At 
the nozzle exit, several distinct regions of the flow can be defined. 

The segment MQ is dependent on the vertical expansions and as well as 
on the lateral expansion. If the lower cowl length (L) is sufficiently long, 
some of the expansion waves generated at A will be reflected and the region 
MO is additionally affected by this reflected wave (where the position of 0 
may be above or below that of Q) . 

The first reflected wave from the upper wall reaches the exit station at 
point H, hence the flow between Q. and H (or 0 and H) undergoes the dual ex- 
pansion cr^ and, ^2 but receives no reflected waves. The waves generated or 
reflected at the upper walls cross the exit line between H and P and have an 
important effect on the vehicle pitching moments. Note that these waves can 
be either expansions or compressions. The value of the turning at the cowl 
(o£) * s a function of and the cowl length L as well as the amount of re- 
compression introduced between B and M. 

Consider the optimum nozzle that can be obtained for a fixed area - i.e., 
the lengths, shapes and turning angles are arbitrary. Optimum conditions 
require that the exit flow be parallel to the free stream and the exit pres- 
sure be constant for a uniform initial flow field. Such a nozzle is one 
where the region OH extends across the entire exit flow station, thus, the 
cowl and vehicle lengths and recompressions are such as to cancel the expan- 
sions emanating from A and B. The turning angles c^ and o 2 must be selected 
such that " c 9 - 3 where 3 is the initial flow inclination with respect 
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to the external flow field. Such an "ideal" nozzle will generally result in 
excessive cowl and possibly vehicle lengths. 

For the real nozzles considered, the expansions produced at A and B wi 1 1 not 
be fully captured thus producing a nonuniform exit flow. In the uniform 

region OH, the expansion results from a turning through + 0 £ plus the 

lateral expansion, and the flow deflection angle equals 3 + ~° 2 ‘ ' n 

regions HP and OM, the resultant expansionis dependent on the percentage of 

the waves from A and B that are reflected and the recompression introduced 
in these regions. The resultant reflected waves can be either expansion waves 

or compression waves. Note that if the reflected waves are compression waves, 

they tend to decrease the exiting velocity and the local Mach numbers, yield- 
ing a maximum Mach number in the region OH. 


For those nozzle flow fields where a uniform region such as OH occurs at the 
exit station, the flow field between M and H is independent of the vehicle 
geometry between A and P. A value of is assumed and for this value, sever- 
al values of a 2 are considered. The f i rst down-runn i ng characteristic from 
A (ART) determines the minimum cowl length required to obtain reflection at 
the lower cowling. This line is first obtained two dimensionally ([0-y]^ = 
( 8 -y)^; [ 0 -y]^T - 63 " °2 " w ^ ere Py represents an expansion by a turning 
angle of - 3. After obtaining the two dimensional position of point T, a 
correction is introduced in terms of the area ratio due to lateral expansion 
between B and T. This correction is introduced as follows: the Mach number 

obtained at T from two dimensional considerations defines an area ratio in 
expansion from the burner Mach number to My. This value is multiplied by the 
value, of the lateral expansion ratio (z^.z.) and this new expansion ratio 
yields a new (larger) Mach number M| and the line RT is redrawn employing the 
value of y| based on the modified Mach number. 


Employing this procedure, several other characteristic lines corresponding to 
different values of turning at point A ranging from 0 turning to 0 ^ may be 
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determined and thus the lines AL, AM and AQ. may be generated, with similar 
procedures used to obtain the corresponding reflected waves. 

From the wave diagram constructed, the value of thrust for : the mass contained 
between M and H may be determined for several values of , a , L and recom- 
pression between T and M. The thrust obtained depends only on these para- 
meters and does not depend on the shape of the vehicle undersurface between 
A and P. The calculation may be first performed for a straight lower cowling 
( BM straight). After, some wave cancellation can be easily introduced for 
each value of L between T and M and the optimum value obtained by trading off 
a reduction in expansion for a decrease in flow angle at the exit station. 

For each value of cr^ chosen, a curve is obtained yielding the variation of 
thrust (Th) as a function of cr^ for a given value of L and, hence, a curve 
of c.j vs. L with thrust as a parameter can be generated. 

In order to optimize the calculation for the values of cr^ , L and a ^ selected 
the following possibilities are investigated; zero reflection and cancella" 
tion, total reflection and half reflection. For each case, the position of 
M is determined and the external drag is obtained. Then, the thrust of the 
mass crossing MH minus the drag is plotted in terms of percent reflection and 
the optimum value obtained for each set of L and a. The loss in thrust in 
percent as compared to the optimum is used as a parameter for comparison. 

Now, the flow between H and P may be analyzed. For this analysis, the values 
of a.j selected for the analysis of the lower region are considered for a given 
value of a.^) hence the variable is the curvature of FP. For each variation of 
the height M, an equivalent variation of the height of P is chosen which then 
defines the curvature of FP (if a parabolic recompression is selected) and the 
thrust of HP is obtained. 

After performing this analysis for several values of c 2 , the range of para- 
meters which require more accurate analysis are defined and approximate values 
of forces are also obtained. This approach is quite simple since changes in 
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o j and o ^ require only the addition of several waves. The same is true for 
changes in reflections. 
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SECTION I ! I 

SOURCE FLOW CHARACTERISTIC APPROACH 

Step (1) of the design process has defined the overall nozzle expansion ratio, 
the transverse to longitudinal expansion ratio and the approximate cowl length. 
This next phase serves to refine this data and fill in the details of the noz- 
zle shape. This approach will be based on quasi two dimensional considerations. 

In the preliminary design of the vehicle undersurface and cowl the quasi two 
dimensional approach must consider, at least to a first approximation, the 
lateral extent of the nozzle. Consider the following: the basic NASA 

scramjet/ramjet design philosophy assumes multiple injector burner configura- 
tions. The burner flow can be separated into regions with combustion occur- 
ring at different longitudinal stations; expansion of the flow begins in the' 
combustor and continues into the nozzle. The resulting wave pattern directly 
influences nozzle performance. 

Figure (6) is a schematic of a multiple burner design which results in a 
lateral distribution which may be roughly approximated by a line source type 
flow, which more closely approximates the physical case than the usual two 
dimensional conditions assumed, which correspond to single injector designs. 

The burner flow to a first approximation (neglecting shock waves) can be 
represented by an average velocity, but gradually changing direction, as in 
a source type of flow. 

In the present approach, the lateral extent of the nozzle will be treated as 
a cylindrical source term in the first pass design of the vehicle undersurface 
and cowl. This produces a three dimensional effect, but does not account for 
side wall curvature and shock losses. However, during the course of the calcu- 
lation it will be possible to specify, at a given axial station, a change in 
the origin of the source flow. This corresponds to a change in the lateral 
v/all angle a as shown in Figure (7). 
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The same flow technique is employed in conjunction with a two dimensional 
second-order characteristic procedure capable of analyzing the aerodynamic 
performance of typical nozzle configurations selected from simplified anal- 
ysis. Such a configuration is depicted in Figure (8). 



FIGURE 8. 'TYPICAL SCRAMJET NOZZLE 
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However, the calculation procedure is not limited to these configurations 
but can be readily adapted to calculate other two dimensional configurations. 
This generality results from the use of three coordinates systems, axisym- 
metric, axially expanding (source type flow) and Cartesian (plane two dimen- 
sional), as described in detail in references (2) and (3). Automatic pro- 
visions for switching from axially expanding to Cartesian coordinates at a 
specified axial station and multiple source origins are provided for. These 
features allow the lateral nozzle area variation to be accounted for in a 
quas i -two td i mens iona 1 manner. A higher order calculation would involve a 
fully three dimensional calculation which would locate the lateral waves. 

The working fluid is assumed to be a hydrogen-air mixture in frozen or 
chemical equilibrium. 

The following boundary conditions are provided for in the calculation. 

(1) Wall boundaries 

(2) Shock boundaries 

(3) Contact surface 

(4) Underexpans ion interaction 

(5) Overexpansion interaction 

(6) Prandt 1 -Meyer 

A running integration of pressure on the nozzle surfaces is performed which 
yields thrust, lift and pitching moment in vehicle coordinates. 
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SECTION IV 

ANALYSIS OF ENGINE-VEHICLE INTERFERENCE EFFECTS 


For our purposes engine-vehicle interference effects will be classified into 
two categories; 

a) Forces due to the presence of the engine plume 

b) External cowl and/or fence forces 

The establishment of suitable engineering criteria for the computation of 
these forces, within the context of a parametric nozzle analysis, are derived 
from linearized supersonic flow-. Conical flow relations are employed in 
corner regions of the external flow, and use is made of the concept of preser- 
vation of lift as an interference effect The latter results in a very 

useful integral relation which is applicable to nonlinear flow fields. The 
conical nature of the external cowl flow field, yields a closed form solution 
for the corner pressure distribution. In addition, a simple correlation for 
the interference pressure in an internal compression corner has been developed. 

A . Plume Interference Analysis Preservation of Lift as an Inter- 
ference Effect - Consider a wing, as discussed in reference (7), with super- 
sonic edges whose trailing edge is perpend icular to the flow direction as in 
Figure (9a). The lift coefficient due to the deflected element has the two 
dimensional value ha (M = *./Y) , when based on the area of the deflected 

O 00 

element. This concept may be viewed as an interference effect, since the 
undeflected part of the wing, while generating no lift, helps preserve the 
lift of the deflected part. Now consider the wing in Figure (9b). If the 
undeflected area to the right of CF were removed, the edge CF would be sub- 
sonic and (for a Q > 0) the lower surface affects the upper surface via an 
upwash field. Thus, the lift in the triangle CEF would only have half of 
its two dimensional value. If the flat plate at zero angle of attack is re- 
inserted it gets a positive apparent angle of attack due to the upwash field, 
and will carry positive lift. Thus the interference is favorable. Further- 
more, the interference is mutual; the lift due to the apparent angle of 
attack on the right half will induce additional lift on the left half, etc. 
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The final result of this mutual interference is that all the lift is preserved 
In using the word "all", the two dimensional value is the standard of com- 
parison. However, only part of the lift is carried by the deflected part. 

The remainder is carried by the undeflected part. 



FIGURE 9 
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Consider the underexpanded nozzle exhaust flow field shown in plan view 



FIGURE 10a 

is the local undisturbed external flow, Region II is the disturbed external 
fl ow. Regions III, and IV are the disturbed exhaust flow and Region V is 
the undisturbed exhaust flow. The normal force is defined as the integral 


E /P.dA 
i n 

i 


= F 


Now assume that the nozzle exhaust were separated from the undisturbed vehicle 
flow by a plate ab as in Figure (10b). The normal force is defined as the 
integral 

E P.dA = F 
I ,V / 1 n n 2 


A. 
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FIGURE 10b 

From the theorem of Preservation of Lift F j = F^. This may be illustrated 
by the following simple example. Define the flow in Region I, and V as 

P v - 10 Pj P , “ 1 ' 

My » 3 M, = 8 

Y v 13 1-25 - Yj “ 1 .4 

A y = A l A ! = 1 

Then ? ^ ~ 11 Pj. Using the above conditions and computing the integral of 
the pressures in Regions i, II, III, IV, and V, yields a value of F j - 10.9 P 
Similar results were obtained by performing the indicated computation for 
other arbitrary flow conditions. 

Thus the concept of preservation of lift yields the following general result: 
"The integrated force for the complete flow field can be determined from the 
equivalent two dimensional pressure distribution; i.e., replacing the plume 
by plate ab and computing the equivalent two dimensional flow yields the same 
integrated force". 
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However, since the plume is three dimensional an additional concept relating 
the integral of the disturbance pressure on the vehicle to the integral of 
the disturbance produced along the vertical plate (see figure 10c) is re- 
quired in order to apply the preservation of lift theorem. That is the con- 
cept must be generalized to disturbances originating on the plate above the 
vehicle undersurface. 



Side View of Plate 
FIGURE 10c 

The extension of these concepts to three dimensional flow fields is discussed 
below within the framework of linearized flow. 


Integrals of source induced flow fields : Consider a semi - i nf i n i te horizontal 

elemental surface QR having width Sx aligned in the streamwise direction as 
in Figure (11). The strip is at a distance £ from the x~y plane and has a 
source intensity f(^) associated with a disturbance pressure AP. The trace 
of the Mach cone from Q. on the x,y plane is denoted by ACB. The intersection 
of the Mach plane from point p with the x-y plane is the line AB. From 
Reference (8) the following integrals are shown to be valid. 


<5udy = 


U t (P)6w 


(iff.(P)/g)6to 


B j 6wdy = W t (P) <$w = - f (P) 6co 

JA 
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The equations thus yield useful insight into the lateral spreading of dis- 
turbances in three dimensional flows. If the distance 5 # 0 then the above 
results apply even if the surface element is rotated about its streamwise 
axis. if in the limit f(£) 610 F(^) then from the above and Figure ( 1 2 ) 

we obtain 

B 
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FIG. 12 Integral relations. 


Replacing the source distribution by the disturbance pressure associated with 
the presence of plate ab leads to the following conclusion. The integral of 
the disturbance pressure on the vehicle is equivalent to the two dimensional 
values; i.e., since the disturbances on either side of the plate are equal 
and of opposite signs the resultant interference effect is equivalent to a. 
quasi two dimensional computation in the presence of plate ab. This is the 
three dimensional counterpart of the simple two dimensional result illustra- 
ted . previously. 

For linearized flow the above result is exact since disturbances on either 
side of the plate propagate at the same free stream Mach number. As the Mach 
numbers can differ substantially on either side of the plate (see Figure 13) 
corrections must be introduc ed for these differences. This is readily 
accomplished since the shape of plate ab is known and the Mach lines are given 
by the characteristics solution. 
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Shock 



Plan View 


FIGURE 13- 


Thus the plume interference problem has been reduced to an equivalent quasi 
two dimensional problem except for the correction introduced above. 

B. Conical Flow Solutions - The cowl can be represented by inter- 
secting planar surfaces forming corners with included angles of 270 ° (external 
corner) or 90 ° (internal corner) as in Figure (1*0 4 The flow field in each 
corner is conical, and as will be demonstrated, classical conical flow 
solutions yield the qualitative nature of the solution, but must be corrected 
for nonlinear effects introduced by local flow conditions. 
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Vehicle 



FIGURE 1*1. COWli AND FENCE 

Classical conical flow is one in which the flow properties are constant along 
rays t = Xj/X^ as in Figure (15)- If the Xj X^ plane defines a wing of 
finite span with deflection <5 and free stream Mach number > 1 , the flow 
on the tip region (shaded area, Figure 15) is given by the conical solution, 



-1 

cos 


/ E-COS Ip \ 

v 1- e cosip' 


If beyond the wing the boundary condition <j> n = o (i.e., disturbance velocity) 
is applied to a plate rotated through an angle i ■ v as in Figure (15b), then 
the solution becomes 



. \ T\ , ' D 

\ -A A 

/\ ' v / \ 
'1C . 

rr\r\ \ 


conical 
reg i on 


Sweep Angie 

-1 

tan t 


sin" 1 (--) 


No Dihedral 


13 0 X \ 
/ ' \ 

tu 


No Dihedral (a) 


$ <=> 6 
n 


cos » tan <p tan \i 
e = 3t 


_ 26 -1 r/ £-cos \p 

C a cos [ 1 

p TT0 1— r rn? ih' J 


' 1 -E COS V 


e /"/V "A 

1 if ■« •» 

Wi th D ihedral (b) 


— oos" 1 [-(f2*jL); 
it 3 1 -e cos ifi 1 


L K <5 0 

n 2 


V, i 


<J> 6, 

n 1 


Wing wi-th Dihedral (c) 


f. f. r -1 ,C"C.OS ill' 

C » L [ir-cos (———I — 

P • P 1 ~e cos ili 

Base 

, , ,, -1 ,e+cos ill' v, 

4- 5-/6. cos ( ■ ■■ — )] 

‘ 1+e cos ip 1 


FIGURE 15 


C ** 
p 

Base 


v = d ihedral angle 
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Cp 


26 

irg 


-1 

cos 


(E^COSJ. ,, 
1- e cos ip 


where ip = (t r/v)^ 


This is a so called dihedral solution, the details of which can be found in 
References (9),(10) and (l l) , if the boundary condition = 0 is replaced by 

$ n = ’"*2 t ^ ien ^ superposition of linear solutions 


Cp = C 


p Base 


„ 

1 1 

-1 

/ e - cos ip s j 

7f “ COS 

\ ^ ) ] 
. 1 - e cos ip J 



'Y^ase* cos 


-i- cos ip^ ^ 
1 + e cos ip^ 


This represents the solution for a wing with dihedral angle v. Figure (16) 
is a plot of the above equation for v - 270° (external corner). These re- 
sults are valid for the parameter M 6< < 1, i.e., linearized supersonic 
flow. For ^ 0(1) the results can be improved if (^^base^ ' s re Pl ac ®d 
by the ratio of the local 2-D pressure coefficients. For the hypersonic 
approximation tangent wedge solutions are appropriate for these values. 
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For the internal corner v = 90 the linearized result reduces to 


P-P 


^ ’ ‘VbaJ 


base 00 


However, Reference (12) indicates that this result considerably under- 
estimates the interference pressure. For equal wedge angles, 6 ^ = < ->| ja r e » 
if the linearized interference pressure is viewed as a perturbation from 

the base pressure P, then 
r base 


P-P 

Base = 1 
^Base 


o 

Which agrees with the result of Reference (12) for M to = 3.17 and 6 = 12.1 . 
From the data of References (12) and ( 13 ) an approximate correlation of the 
interference pressure for an internal compression corner has been determined, 
i . e. : 



For M <$, >>1 using the local Mach number rather than M- correlates the 

CO Q 00 

results well even up to M = 20. 

Given the above local corner solutions the cowl forces can be computed by 
integrating over the appropriate cowl areas defined by the plane of symmetry 
and the local Mach cones in the corner. 
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SECTION V 
CONCLUSIONS 


The current ATL parametric approach to scramjet nozzle design allows the 
propulsion specialist to rapidly and logically assess the effects of various 
parameters on engine and vehicle performance. The three step method of in- 
creasing complexity permits one to examine only the most promising configura- 
tions selected from the previous step. Thus extraneous and time consuming 
calculations are eliminated. As part of the above parametric approach ATL 
has developed, in References (l~5), and has delivered to NASA computer programs 
based on the level of sophistication required in each step of the optimization 
process . 

The concept of preservation of lift and the use of equivalent two dimensional 
solutions allows for rapid assessment of the external plume interference 
problem. In addition, the conical flow solutions developed, significantly 
simplify the force computation associated with external cowl and fence sur- 
faces . 

The integrated parametric methodology described herein 
leading to the optimum vehicle-propulsion system based 
straints imposed by the designer. 


is a logical process 
on the physical con- 
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CURRENT PROGRAM CAPABILITIES 


WAVE PROGRAM 


SOURCE FLOW CHARACTERISTICS 


THREE DIMENSIONAL PROGRAM 


Advantages 


(1 ) Ease of use 


(2) Rapid Calcu- 
lat ions 


(3) Mechanized 


Disadvantages 


(1 ) L imi ted to un i - 

form. initial 
prof ! 1 e 

(2) ! sentropic, ho- 
mentropic, 
constant ® 

(3) limited to 
certain wave 
conf igurat i ons 


(4) Small storage (A) No external 

flow interac- 
tion 


Advantages 


(5) Equ i 1 i br i urn 
H^-Ai r or 

ideal gas 


(5) Quasi 2-D 



D i sadvantages 


to use 

(2) Complex Logic 


(3) Longer running 
t ime 


(k) Not mechanized 
ry (one case at 
a time) 

(5) More storage 


(6) quasi 2D 
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